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Abstract 


A three-dimensional model has been developed for the free-breathing direct methanol fuel cell (DMFC) assuming steady-state isothermal and 
single-phase conditions. Especially the MeOH crossover phenomenon is investigated and the model validations are done using previous cathodic 
current distribution measurements. A free convection of air is modelled in the cathode channels and diffusion and convection of liquid (anode) 
and gaseous species (cathode) in the porous transport layers. The MeOH flow in the membrane is described with diffusion and protonic drag. The 
parameter y in the model describes the MeOH oxidation rate at the cathode and it is fitted according to the measured current distributions. The 
model describes the behaviour of the free-breathing DMFC, when different operating parameters such as cell temperature, MeOH concentration 
and flow rate are varied in a wide range. The model also predicts the existence of the experimentally observed electrolytic domains, i.e. local 
regions of negative current densities. Altogether, the developed model is in reasonable agreement with both the measured current distributions and 
polarization curves. The spatial information gained of mass transfer phenomena inside the DMFC is valuable for the optimization of the DMFC 


operating parameters. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Conventional batteries are soon becoming inadequate for the 
increasing power requirements of portable electronic devices 
and computers [1]. The alternative power supply in the future 
could be the direct methanol fuel cell (DMFC). The DMFC is 
a polymer electrolyte fuel cell (PEFC), where the anode reac- 
tant is a liquid mixture of MeOH and water. The use of MeOH 
gives some advantages compared to the hydrogen systems: the 
storage of the fuel is easier and the refueling systems can have 
a simpler structure. Because of simplicity, space saving and 
economical reasons, some DMFC applications will probably 
be totally passive-feed DMFCs, free-breathing and without any 
external pumps. 
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In order to optimize performance values for the DMFC, it 
is important to have even current and temperature distributions 
inside the fuel cell. This can be achieved, when the concen- 
trations of the reactants have optimum values. In many cases 
the goal of the fuel cell development is also to achieve an even 
current density profile under a variety of operating conditions. 
As a consequence, current distribution measurements are very 
valuable for the better understanding of the complex relation- 
ships between different variables and operating parameters of 
the DMFC. The current distribution measurements have been 
used in many PEFC [2-11] and DMFC [12-16] studies during 
recent years. 

In the case of the DMFC, the MeOH crossover (MeOH per- 
meability through the membrane) is a very dominating loss 
mechanism [17] and it has been investigated widely [18-29]. 
However, most of these studies have been performed with- 
out current distribution measurements. The MeOH permeating 
from anode to cathode poisons the Pt catalyst and causes mixed 
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Nomenclature 


concentration (mol m7? or M) 

diffusion coefficient (m? s7!) 

voltage (V) 

Faraday’s constant (96 485 C mol!) 
acceleration due to gravity (9.81 ms~7) 
current density (A m~?) 

thickness (m) 

molar mass (kg mol`!) 

molar flux (mol m~? s7!) 

pressure (Pa or atm) 

universal gas constant (8.314 J mol`! K7!) 
relative humidity (%) 

temperature (K or °C) 

velocity (m s7!) 

diffusion volume 

flow rate (m? s7!) 

horizontal coordinate in the model (m) 
mole fraction 

vertical coordinate in the model (m) 
coordinate perpendicular to xy-plane (m) 


Subscripts and superscripts 


reference 

anode 

ambient 
average 

cathode 

fuel cell 

critical 

effective 
electrode 
species 

species 

liquid 

protonic phase 
membrane 
porous transport layer 
reference 
electronic phase 


Greek letters 


transfer coefficient of electrochemical reaction 
porosity 

overpotential (V) 

intrinsic permeability (m7) 

electro-osmotic drag coefficient 

dynamic viscosity (kg m7! s7!) 

density (kg m~?) 

conductivity (S m7!) 

potential (V) 

parameter for MeOH oxidation rate at cathode 
(ms™') 


potentials, which decrease the fuel cell performance [30,31]. 
The polarization losses are caused mainly because MeOH is 
oxidized heterogeneously in the presence of oxygen [32]. Ye 
and Zhao [33] have observed a significant decline of open cir- 
cuit voltage (OCV) in the DMFC, when the oxygen flow rate is 
reduced below a critical value. This can be explained with the 
coexistence of galvanic and electrolytic reactions in the DMFC, 
even when the cell is operated under open circuit conditions 
[34,35]. These reactions produce self-discharging currents in the 
DMFC causing MeOH consumption and H? evolution. On the 
contrary, when the oxygen flow rate is above a certain value, the 
OCV is observed to be very insensitive to this quantity [33]. For 
investigation of that kind of phenomena, the DMFC with a seg- 
mented cathode can be useful and provide detailed information 
on local concentrations of the reacting species, on which part of 
the DMFC most of the current is produced and how much the 
MeOH crossover phenomenon decreases the cell performance 
locally. 

Many fuel cell models for the PEFC and the DMFC have 
been developed, see, e.g. review articles of Wang [36], Weber 
et al. [37] and Yao et al. [38], but most of the DMFC model 
validations have only been done using measured polarization 
curves, e.g. [38—44]. Models for passive-feed DMFCs with heat 
transfer effect have been set up by Chen and Zhao [45](1D) and 
Rice and Faghri [41](2D). Wang and Wang [43] have proposed 
a 2D two-phase model and Danilov et al. [44] a 3D two-phase 
model of a liquid-feed DMFC. Kjeang et al. [46] have published 
a 3D model of a flowing electrolyte DMFC to simulate MeOH 
crossover by convection and diffusion and Liu and Wang [47] 
a 3D two-phase DMFC model to investigate electron transport. 
Ge et al. [48] have made a 3D single-phase model of the DMFC 
to study the effects of MeOH crossover, porosities, MeOH flow 
rates and channel shoulder widths. 

Ju et al. [49] have set up a 3D PEFC model (H2 feed) and val- 
idated it with measured current distribution data. They also con- 
cluded that model validations against only polarization curves 
is insufficient and often misleading [49]. Because the chemical 
reactions and connections between different operational param- 
eters in many cases are more complicated in the DMFC than in 
the hydrogen fed PEFC, there is a true need for DMFC model 
validations using current distribution data. However, this kind of 
3D model for the free-breathing DMFC has not been reported so 
far. In this paper a 3D model for the free-breathing DMFC is pre- 
sented and the model validations are done using our previously 
published current distribution measurements [12]. 


2. Modelling 
2.1. General model description, geometry and assumptions 


A three-dimensional model of the free-breathing direct 
methanol fuel cell (DMFC) has been developed to investigate 
mass transfer phenomena. In particular the MeOH crossover 
phenomenon is investigated and model validations are done 
using the current distribution measurements [12]. The model 
is implemented using the commercial Comsol Multiphysics 3.3 
program (Comsol AB). 
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In this model, the Navier-Stokes equation (Eq. (1)) is used 
for modelling Newtonian flow of MeOH and water in the anode 
flow field and also free convection of air in the cathode flow 
field. The Stefan—Maxwell equation (Eq. (8)) is used for mod- 
elling multicomponent diffusion of gaseous species in the porous 
cathode transport layer (PTL) and Brinkman’s equation (Eq. (2)) 
for modelling the flow of MeOH and water in the anode PTL. 
The Tafel equation (Eq. (17)) is used to describe electrochemical 
reactions at the cathode and the linearized Butler—Volmer equa- 
tion (Eq. (18)) at the anode. The MeOH flow in the membrane 
(MeOH crossover) is described with Fick’s law of diffusion and 
drag of protonic current (Eq. (19)). The only fitted parameter 
in the model is the parameter y for the MeOH oxidation rate 
at the cathode. y is fitted according to the measured current 
distributions. 

The model geometry is illustrated in Fig. | and the main 
flow directions of MeOH-—water solution and free air convec- 
tion are also marked in the figure. The air channels of the cell 
(y-axis) are vertically oriented. Producing the DMFC model 
in 3D causes some restrictions and compromises to the mod- 
elling work, because the model solution should be found in a 
reasonable time with the available computer resources. Much 
more detailed models can be done in 1D [40,45], 2D [41,43] 
or if only some part of the DMFC, e.g. the anode is modelled 
[50,51]. In our case, the available computer resources practically 
means that the number of elements and degrees of freedom in 
the model should not exceed 10°. The main simplifications to 
the real DMFC geometry, which is described in detail in [12], 
are that the electrodes are considered only as boundary condi- 
tions and the whole area of the PTL (z-plane) is considered as 
a current collector instead of the segmented electrode used in 
the experimental work. Also, the flows of reactants (both the 
inflows and outflows) are set to be only in the main flow direc- 
tions (see arrows in Fig. 1) instead of in the separate air channels 
of the cathode and the serpentine parallel MeOH channels in the 
anode of the experimental DMFC. The geometry can also be 
scaled down (1:10) in the x-direction by adjusting the corre- 
sponding flow parameters, which leads to a significant decrease 
in computing time without affecting the modelling results. The 
following assumptions are used in the model: 


MeOH flow 


Air convection 


Fig. 1. The modelled geometry and the main flow directions of reactants. 


ee 


. Steady-state isothermal and single-phase conditions. 

2. No phase transitions (gas condensation or liquid vapouriza- 
tion) between liquid water (anode side and membrane) and 
water vapour (cathode side). 

3. Ideal gas assumption for gaseous species and all flows are 
assumed to be laminar. 

4. Liquid species (MeOH and water) are considered to be New- 
tonian fluids. 

5. Complete MeOH oxidation (no partial oxidation or side reac- 
tions). 

6. No gas fluxes through the membrane. 

. CO2 and H; fluxes are neglected. 

8. All electrochemical reactions take place at indefinitely thin 

electrodes and no reactions occur in the membrane, porous 

transport layers or flow field regions. 


N 


2.2. Governing equations 


2.2.1. MeOH flow field and porous anode transport layer 
The flow of MeOH and water at the anode flow field can be 
described with the Navier-Stokes equation: 


du 


v 1 
B= Pe, (1) 


pu- Vu = uV°u 


where jz is the dynamic viscosity of the MeOH—water solu- 
tion. The flow through the porous anode transport layer (PTL) 
can be described according to Brinkman’s equation (Brinkman’s 
extension of Darcy’s law): 


uV?u— Zu =Vp (2) 


where x is the intrinsic permeability of the porous medium. At 
the anode flow field and the PTL, the molar fluxes of MeOH 
and water are continuous (V - N; = 0) and are modelled with 
the equation for convective diffusion: 


N; = —D;Vc; + cju (3) 


2.2.2. Air flow field and porous cathode transport layer 

The developed three-dimensional model for the air flow field 
and the cathode PTL is based on a 2D model developed by Men- 
nola et al. [52]. The steady-state continuity equation for mass 
(V - pu = 0) is used and the free convection of air at the cath- 
ode channels is described with the Navier-Stokes equation (Eq. 
(1)). The buoyancy in the vertical cathode channels (y-direction) 
is caused by density differences of air resulting from thermal 
expansion (Toe > Tamb) and changes in gas composition (elec- 
trochemical reactions). The pressure variations in the flow field 
region are assumed to be small, i.e. Eq. (1) simplifies to the form: 


uV’u +g — po) =0 (4) 


where the density of air p is calculated according to Eq. (5): 


p =c XM; 6) 
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The total concentration c can be calculated from the ideal gas 
law: 


ae ee 
RT 


The molar fluxes of chemical species are continuous (V - 
N; = 0) in the air flow field and are defined according to Eq. 
(3). The mole fraction of nitrogen is calculated from the mole 
fractions of oxygen and water vapour: 


c 


(6) 


Xn, = 1 — Xo, — Xm,o (7) 


In the porous cathode transport layer (PTL) the 
Stefan—Maxwell equation is used for modelling multi- 
component diffusion: 


Wy = 


~ cDij 
ps 


XN; — XjNi) (8) 


where Dj; is the binary diffusion coefficient of the gas pair i — j. 
Nitrogen is inert and its flux is set to zero. The molar fluxes for 
oxygen and water vapour can be solved according to Eq. (8): 


Do 3,820 PN>,0» 


No, = —c 
DyN, 0, Xm o + Dni, m oXo, + Dor,H20XN> 
Xo, DN3,H20 
(FES Y yio + VXo,) + VXor (9) 
XN, Do,,H,0 
Do),4 0 PNy,H0 
Ny,0 = 


C x 
DN, ,0.XH,0 + DNz,H,0X0, + Dor,Hy0XN> 


(Serta 


(VXm o + VXo,) + VXino } 
Xn, Do,,H,0 ? . ° 


(10) 


2.2.3. Electrode reactions 

According to our experimental results [12] and observations 
of, e.g. Ye et al. [34] and Kulikovsky et al. [35], two differ- 
ent regions depending on local conditions inside DMFC can be 
distinguished: 


1. Normal DMFC operation (Galvanic cell): 


anode: CH30H + H20 —> CO) + 6H* + 6e” (11) 


cathode: 302 + 6H* + 6e7 > 3H20 (12) 
2. Electrolytic cell: 


anode: 6H* +6e~ — 3H» (13) 


cathode: CH3;0H+H 20 > CO,+6H'+6e (14) 

Reaction (12) is the main reaction at the cathode of the 
DMFC, but because of the MeOH crossover phenomenon, 
MeOH oxidation will also occur at the cathode and the cor- 
responding anode reaction is then hydrogen formation. When 
the MeOH crossover rate through the membrane is low, reaction 
(12) is dominating at the cathode. In this model, the cathode 


and anode electrodes are taken into account as boundary condi- 
tions and electrochemical reactions take place at indefinitely thin 
electrodes. The drain and source terms of MeOH and water are 
included in the model according to the anode and cathode reac- 
tions, but the fluxes of CO7 and H3 are neglected. In the DMFC, 
the cathode and anode PTLs act also as current collectors for 
electricity and only an electronic potential exists: 


V- (-08;V¢s) = 0 (15) 


The transfer of electronic current into protonic current takes 
place at the boundaries between the anode and cathode PTLs and 
the membrane. The current production at the cathode (reaction 
(12)) is approximated with the Tafel equation: 


f f &cF 
lee = iĝ exp (- a ne) (16) 


where the cathodic overpotential ne = ¢; — ¢m — E°. Because 
the MeOH crossover phenomenon decreases the total current 
production at the cathode, the following modifications to Eq. 
(16) are introduced: 


. < [co Qc F 
fee = iĝ (2) exp (- F ne) — 6FcMeon Y (17) 


O2 


where the parameter y describes the MeOH oxidation rate at the 
cathode. This parameter is to be fitted according to the measured 
current distributions. The effect of local oxygen concentration is 
also taken into account in Eq. (17) and the dependency between 
MeOH oxidation rate and MeOH concentration is assumed to 
be linear. When the DMFC is operating in galvanic mode, with- 
out any MeOH crossover, the Tafel equation (Eq. (16)) is the 
most reasonable way to describe current production at the anode. 
The drawback of the exponential Tafel equation is that it cannot 
be used alone without any additional terms or fitting param- 
eters, in order to describe also the local electrolytic domains 
found experimentally [12]. Thus, a major simplification for the 
anodic current production is made by assuming only low anodic 
overpotentials and approximating the anodic current production 
(reaction (11)) with the linearized Butler-Volmer equation: 


e n | CMeOH \ QaF 
a= a/ L ) T Na (18) 


where the anodic overpotential na = m — ¢s. The effect of local 
MeOH concentration is also taken into account in Eq. (18). 


2.2.4. Membrane 

Noreactions occur in the membrane and only a protonic phase 
exists (V - (—o;..,V¢m) = 0). The MeOH crossover through 
the membrane occurs because of electro-osmosis and the gradi- 


ents of pressure and concentration: 


AMeOH. K 
MeOH = = i+ p MeOH Y P — DMO VC (19) 


where Ameou is the electro-osmotic drag coefficient (number of 
MeOH molecules dragged by each proton). The flux induced by 
the pressure gradient is assumed to be small compared to other 
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fluxes and so it is neglected in the model. The last term in Eq. (19) 
is the diffusion flux according to Fick’s law of diffusion. Because 
of relatively high water content on both sides of the membrane 
(liquid water at the anode and water vapour and water production 
at the cathode), the membrane is assumed fully hydrated and 
water flow in the membrane is modelled only with the drag of 
protonic current: 


AO . 
Ngo = Ea (20) 


2.3. Boundary conditions 


2.3.1. MeOH and water 

The concentrations and velocities of MeOH and water are 
fixed at the beginning of the anode flow channel. Whereas at 
the end of the anode flow channel, the concentration gradients 
are set to zero and the pressure is fixed to the constant pres- 
sure of ambient air. The boundary conditions between the anode 
PTL and the membrane are set according to MeOH and water 
consumption (Eqs. (11) and (18)) and the boundary conditions 
between the membrane and the cathode PTL according to Eqs. 
(14) and (17) (oxidation of permeated MeOH). 


2.3.2. Oxygen and water vapour 

The input concentrations in ambient air are calculated accord- 
ing to Eq. (6) and the vertical velocity gradients are assumed 
to be zero (fully developed flow) at the end of the air chan- 
nel. The boundary conditions between the cathode PTL and the 
membrane are set according to oxygen consumption and water 
vapour production (Eqs. (12) and (17)). Gas fluxes to the mem- 
brane are set to zero and phase transitions between liquid water 
and water vapour are not taken into account. The pressure is set 
to be constant at both ends of the air channels, because they are 
open to the ambient air. 


2.3.3. Electronic and protonic potentials 

The electrodes are described as boundary conditions in the 
model. The electronic potentials are fixed at the boundaries 
between the flow fields and the PTLs. Because the whole area of 
the anode and cathode PTL acts as a current collector for elec- 
tricity in the model, the electronic potential is set to zero at the 
boundary between the MeOH flow field and the anode PTL and 
to the cell voltage at the boundary between the cathode PTL and 
the air flow field. The protonic potential exists only in the mem- 
brane and the fluxes of electrons and protons are fixed between 
the membrane and the PTLs according to their consumption and 
production (Eqs. (17) and (18)). 


2.4. Parameters and variables 


All parameters in the model are assumed to be independent of 
concentrations of reacting species, e.g. the temperature depen- 
dent dynamic viscosity value for the MeOH—water mixture is 
calculated as an average value of the MeOH concentrations 
(0.5—10 M) used. All common parameters and variables are pre- 
sented in Table 1. The model variables are given the same values 


Table 1 

Parameters and variables 

Parameter Value 

Teel CC) 30, 50, 70 
CMeou (mol dm~?) 0.5, 1,3, 5, 10 
Vmeon mo (m? s7!) 3.3 x 107°, 17 x 107°, 130 x 107° 
lpu (m) 300 x 1076 
lmem (m) 200 x 1076 
My, (kg mol!) 18 x 1073 
Mwcou (kg mol!) 32 x 1073 
My, (kg mol7!) 28 x 1073 
Mo, (kg mol7') 32 x 107° 
Eo (V) 0.6 

Tamb (°C) 25 

RH 50% 

Po (Pa) 1.01 x 105 
på" (kg m~’) 1.177 

Po, (Pa) 20945 

Py, (Pa) 78794 
Pio (Pa) 1585 

ct (mol m~?) 0.136 

re i (mol m~?) 100 

of, (Sm!) 10° 

Kpu (m°) 1x 107"! 
AMeOH = AHO 2.5 

€mem 0.28 

Eptl 0.6 

Ola 0.24 

Oe 0.875 

is (A m~?) 1 

i, (Am7?) 103 

Pcr, (Pa) 3.393 x 10° 
Pcro, (Pa) 5.04 x 10° 
Per,H2O (Pa) 22.1 x 10° 
TerN> (K) 126.2 
Ter,O2 (K) 154.6 

Ter, H20 (K) 647.1 


as they had in the experimental measurements [12]. The den- 
sity of air and partial pressure values of oxygen, nitrogen and 
water vapour are calculated assuming 50% relative humidity of 
the ambient air. The thickness of the compressed porous trans- 
port layer is assumed to be 300 um and the thickness of the 
compressed and fully hydrated Nafion® 117 membrane 200 um. 

The values for the exchange current densities iĝ and iĝ and 
the open circuit voltage Eg are assumed to be independent of the 
cell temperature. The same assumption is also made for the fol- 
lowing parameters: the conductivity of the porous transport layer 
nul is assumed according to Refs. [47,53,54]. The intrinsic per- 
meability «py and the transfer coefficients of the electrochemical 
reactions œa and a, are assumed to be the same as in Ref. [43], 
water drag coefficients AMeou and Ay,0 as they are in Ref. [26] 
and membrane porosity €mem as in Ref. [48]. The value of the 
porosity of the transport layer €py is as the one used in, e.g. Refs. 
[45,48,52,53]. Values of the critical pressure and temperature 
(Per and Ter) are taken from Ref. [55]. 

All temperature dependent parameters are presented in 
Table 2. The conductivity of the membrane is calculated accord- 
ing to Ref. [56], the diffusion coefficient of MeOH in water 
according to Ref. [57] and diffusion coefficient of MeOH in 
the membrane according to Ref. [48]. Before inserting DIKOH 
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Table 2 

Temperature dependent parameters 

Parameter Teen = 30°C Teen = 50°C Teel) = 70°C 
Hair (Pas!) 1.86 x 1075 1.95 x 1075 2.03 x 1075 
umeon (Pas!) 7.7 x 1074 5.3 x 10-4 4.2 x 1074 
Omem (Sm7!) 7.8 10.2 12.8 

Dg, (m? s7!) 2.05 x 1075 2.30 x 1075 2.55 x 10-5 
D$ o (m? s7!) 2.58 x 1075 2.88 x 1075 3.20 x 1075 
DRO (ms!) 1.9 x 107° 3.1 x 10-° 4.7 x 10~° 
prem | (m? s7!) 6.3 x 10710 1.0 x 107° 1.7 x 10~° 
Do ,Ny (m? s7!) 2.12 x 1075 2.38 x 1075 2.66 x 1075 
Do, mo (m? s~!) 2.60 x 1075 3.02 x 1075 3.47 x 1075 
Dn, mo (m? s7!) 2.72 x 1075 3.16 x 1075 3.63 x 1075 


values into Eq. (19), the effective diffusion coefficients are cal- 
culated using the Bruggeman correlation (Eq. (21)) to take into 
account the porosity of the membrane: 


Degg = (€!°)D (21) 


The diffusion coefficients of oxygen and water vapour in air 
are calculated according to Eq. (22), which is based on the work 
of Fuller et al. [58]: 


T3((1/Mj) + (1/M;))!/? 


(22) 
PIED! + Ev) 


Dij = 0.0101 


where the value of the diffusion volume Xv is 16.3 for oxygen, 
19.7 for air and 13.1 for water vapour. The binary diffusion 
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coefficients in the cathode PTL are calculated according the 
Slattery-Bird correlation [55]: 


b 
A(T/ \/Tex,iTer,j) (Pori Per,j) Tori Ter.) (A/M) + A/M 
p 


Dij = 
(23) 


where Ter; and Tor,j are the critical temperatures (K), Der 
and Per,j the critical pressures (atm) and M; and M; molar 
masses (g mol~!) of the species i and j, respectively. The val- 
ues of the empirical constants a and b depend on gas pair 
and for non-polar gas pairs a = 2.745 x 1074 and b = 1.823, 
whereas for a binary diffusion of water vapour with non-polar 
gas a = 3.640 x 1074 and b = 2.334. Before inserting binary 
diffusion coefficients into Eqs. (9) and (10), the effective diffu- 
sion coefficients are calculated using the Bruggeman correlation 


(Eq. 21). 
3. Results and discussion 
3.1. Experimental validation 


The current distribution measurement system for a free- 
breathing DMFC is described in detail in [12]. The cathode flow 
field plate is divided into 48 segments and the current distribu- 
tion is measured with a resistor network. A similar approach 
has been used in PEFC studies by Noponen et al. [5,6]. Most of 
the measurement data used in this modelling work is presented 


Measured current distribution i -2 
—ao—oo“—“rrewx Current density (Am?) Max: 78.367 
78 
77.5 
77 
76.5 
76 
i U 75.5 
o = 
2 ee 75 
Current density 40 i Min: 74.996 
Current density (Am?) Max: 53.024 Current density (Am?) Max:-35.946 
53 -36 
52.5 -37 
52 -38 
51.5 -39 
y 51 y -40 
50.5 -41 
50 -42 
49.5 -43 
49 -44 
x x Ga: 
Min: 48.639 Min: -44.464 


Fig. 2. The measured and three modelled current density distributions with different y values at 70° C (5M MeOH, flow rate 17 x 107° m3 s~!). 
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in [12] and measured using the commercial N afion® 117 mem- 
brane and hot-pressed gas diffusion electrodes for the DMFC 
(E-TEK). 

The current distribution measurements are done using 0.5, 1, 
3,5 and 10 M MeOH solutions and the MeOH flow rates are also 
varied in a wide range from 3.3 x 107° to 130 x 107°? m? s7! 
(0.2—8 ml min” l ). Thus, the MeOH feed stoichiometry is varied 
from 3 to 2500. Because the temperature used and the cell current 
affects the mass transport and electrode kinetics, the measure- 
ments are done at three different cell temperatures (30, 50 and 
70°C) and using different average current densities. 

The parameter w in the model (Eq. (17)) describes the MeOH 
oxidation rate at the cathode and it is fitted according to the 
measured current distributions. Typical values for w are rela- 
tively small, in most of the cases y = 107°? to 1078 ms7!. If 


Current density (A mr?) Max: 127.399 


E, „= 0.45 V,30°C 


x 
Min: 53.922 


Max: 101.333 


100 


Current density (A m?) 


Min: 73.778 


Current density (A m2) Max: 79.257 


IE, ,= 0.45 V,70°C 


71 
Min: 70.86 


x 


too high a value for y is used, the shape of the current dis- 
tribution turns impractical and the current production of the 
cell decreases drastically. An example of this situation is pre- 
sented in Fig. 2, where both one measured and three modelled 
current density distributions at 70° C are presented. The flow 
rate of 5M MeOH is 17 x 107° m? s~!and the w values used 
are 107°, 1078 and 5 x 1078 ms~!. When the value of w is 
107° ms~!, the shape of the modelled current distribution is 
quite similar to the measured one: there are no significant vari- 
ations in the x-direction and the current density decreases in the 
y-direction. It can also be seen that the measured current den- 
sity decreases much faster in the y-direction than the modelled 
one with y = 107° ms~! and with increasing y values, the cur- 
rent density distribution turns to negative values over the whole 
electrode. 
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Fig. 3. The modelled current density distributions at 30, 50 and 70°C. Ece is 0.45 V (left) and 0.5 V (right). The flow rate of 0.5 M MeOH is 17 x 107? m? s7!. 
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Fig. 4. The modelled oxygen concentration at the cathode (above), MeOH con- 
centration at the anode (center) and current density distribution (below) at 70° 
C (1M MeOH, flow rate 17 x 107° m? s~!, Ecen = 0.4 V). 
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3.2. Concentration profiles of the reactants 


During the modelling work, the following general observa- 
tions were done concerning the concentration profiles of the 
reactants. 


3.2.1. Cathode side 

The vertical buoyancy of air is caused by density differences 
resulting from thermal expansion (Tee > Tamb) and increases 
as a function of cell temperature. The highest velocity of air 
is naturally located in the center of the air flow field in the x- 
direction. There are also concentration gradients of oxygen and 
water vapour towards and away from the cathode according to 
reaction (12). As a result, changes in gas composition can be 
seen: the oxygen concentration decreases and the water vapour 
concentration increases as a function of cell current. 


3.2.2. Anode side 

The concentrations of MeOH and liquid water decrease as 
a function of cell current according to reaction (11) and there 
is a MeOH concentration profile along the anode flow channel. 
The shape of this profile depends on the MeOH input concen- 
tration and flow rate. Because of diffusion and electro-osmotic 
drag of protons, there are also concentration gradients of MeOH 
and water through the membrane and MeOH is oxidized at the 
cathode according to reaction (14). If the DMFC is operated 
with a relatively high constant current and with a low MeOH 
concentration and flow rate, most of the MeOH is consumed at 
the anode. As a consequence, there is only a small amount of 
MeOH, which can permeate through the membrane and react 
at the cathode. On the contrary, if the DMFC is operated with 
a low cell current and with a high MeOH concentration, the 
MeOH crossover rate is high. 


3.3. The effect of cell temperature 


The performance values of the DMFC depend quite much on 
the cell temperature and power densities increase as a function of 
temperature, mainly as a result of the increased catalytic activity 
(decreased activation overpotential). On the other hand, some 
part of this gain is lost as the MeOH crossover also increases as 
a function of temperature [27,29]. 

When concentration profiles of reacting species at both elec- 
trodes are compared to the corresponding current distribution, 
it is observed that the oxygen mass transfer limitations decrease 
the cell performance most in the free-breathing DMFC, when 
the cell temperature is 30°C. The buoyancy of air increases, 
when the temperature is increased and this can be seen in the 
concentration profiles of oxygen at the cathode and in the cur- 
rent density distributions. An example is shown in Fig. 3, where 
six modelled current density distributions at three different cell 
temperatures (30, 50 and 70°C) are presented. 

It can be seen from Fig. 3 that mass transfer of air is limit- 
ing the cell performance mainly at low temperatures and when 
the cell voltage is near the OCV. Most of the current is pro- 
duced at the beginning of the air flow field and the current 
density decreases in the air flow direction (y-axis), when the cell 
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Fig. 5. The measured (left) and modelled (right) current density distribution at 70° C (3 M MeOH, flow rate 3.3 x 107°? m3 s~!, Ecen = 0.5 V). 


temperature is 30 °C. At higher temperatures, most of the current 
is produced near the MeOH inlet and the shape of the current 
density distribution turns gradually from the air flow direction 
(y-axis) towards the MeOH flow direction (x-axis) as a function 
of temperature. At the same time, there can be seen a dras- 
tic decrease in the local variations of the current density at the 
electrode. 


3.4. The effects of MeOH and oxygen concentration 


High mass transfer overpotentials can exist in a free-breathing 
DMFC, even at relatively low current densities. This results 
mainly from an inadequate air flow towards the cathode. Dur- 
ing the previous experimental work [12] and this modelling 
work, it has been observed that the local oxygen and MeOH 
concentrations have a significant effect on the current density 
distribution in the DMFC. The MeOH mass transfer is also 
limiting the cell performance especially at high current densi- 
ties, but oxygen mass transfer limitations are highly dominating 
especially in the case of the free-breathing DMFC. Similar 
observations have been done, e.g. by Geiger et al. [16], who 
have made segmented DMFC measurements and concluded that 
the current distribution is highly dependent on the O2 con- 
centration. 

The MeOH crossover decreases the performance of the 
DMFC, especially when using high MeOH concentrations 
(5-10 M), a low MeOH flow rate (3.3 x 107° m*s~!) and if 
the cell is operated at low current densities. When the shapes of 
MeOH and oxygen concentration profiles are compared to the 
corresponding cell current distribution, the limiting mass trans- 
fer process can be indentified. In many cases, it can be seen 
whether the oxygen or MeOH mass transport is the more lim- 
iting one, but sometimes both of them substantially affect the 
shape of the current distribution. An example of this situation is 
presented in Fig. 4, where modelling results for the oxygen con- 
centration at the cathode, the MeOH concentration at the anode 
and the corresponding current density distribution at 70° C are 
shown. The flow rate of 1M MeOH is 17 x 107° m? s~!and 
Ecen = 0.4 V. The mass transfer limitations of MeOH and oxy- 
gen can be seen in different directions: the oxygen concentration 


decreases and limits the cell performance in the y-direction and 
MeOH concentration in the x-direction. 

During the experimental measurements, it has been observed 
that the current densities can vary significantly inside the DMFC 
and in some cases even negative currents have been measured 
[12]. An example of measured local negative currents is pre- 
sented at the left in Fig. 5. The cell voltage is 0.5 V and the flow 
rate of 3M MeOH is 3.3 x 107°? m? s7}. The current densities 
are negative near the MeOH inlet, because the DMFC is oper- 
ated with relatively low current density and the MeOH crossover 
is dominating at those cathode segments. The corresponding 
modelled current density distribution (y = 1.15 x 10-? ms~!) 
is presented at the right in Fig. 5 and it is quite similar to 
the measured one. The model predicts that the current produc- 
tion of the DMFC can be locally negative, when the MeOH 
crossover rate is high near the MeOH inlet. The current den- 
sities increase along the MeOH flow direction (x-axis), which 
can be explained with a decrease of the MeOH concentration 
at the anode. The MeOH concentration decreases due to the 
electrochemical reactions and as a result, the effect of MeOH 
crossover is reduced and the current production again turns 
positive. 


3.5. Model agreement with polarization curves 


After model validations according to the measured current 
distributions, the polarization curves can be obtained by chang- 
ing the cell voltage and calculating the corresponding current in 
the z-direction. An example of calculated and measured polar- 
ization curves and power densities at 70° C is presented in 
Fig. 6. The galvanostatic measurement is made using current 
steps of 1 mA cm? and the flow rate of 1 M MeOH solution is 
17 x 107? m? s™!. 

It can be seen from Fig. 6 that ohmic losses are lower in the 
model calculation compared to the real DMFC system. This is 
quite reasonable, because the whole area of the PTL is used 
as a current collector in the model and the contact resistances 
between the current collector and the PTL are neglected. When 
comparing the measured and calculated polarization curves at 
high current densities, it can be seen that the current production 
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Fig. 6. Measured and calculated polarization curves and power densities at 70° 
C (1M MeOH, flow rate 17 x 107° m3 s™!. 


predicted by the model decreases faster than the measured one, 
evidently due to overestimation of the mass transfer overpo- 
tentials by the model. Also other simplifications of the model 
probably cause that there is not a perfect fit between the mea- 
sured and calculated polarization curves, even though the model 
predicts power densities with reasonable accuracy. For future 
work and to achieve an even better fit, it is recommended that 
the model validations are done both with the polarization curves 
and the current distribution measurements. 


4. Conclusions 


A three-dimensional model is developed for a free-breathing 
direct methanol fuel cell. The MeOH crossover and other 
mass transfer phenomena are investigated computationally and 
the model validations are done using the current distribution 
measurements. Even if the model presented contains several sim- 
plifying assumptions, it describes mass transfer phenomena and 
current production in a free-breathing DMFC, when different 
operating parameters such as cell temperature, MeOH concen- 
tration and reactant flow rates are varied on a large scale. The 
model also describes with reasonable accuracy the existence 
of the measured electrolytic domains, i.e. regions of negative 
current densities. 

The current distribution is observed to be closely connected to 
the local concentrations and mass flows of the reacting species. 
The local power production is observed to be proportional to 
the local oxygen concentration at the cathode side. As a con- 
sequence, inadequate air convection together with the MeOH 
crossover phenomenon decreases the cell performance. At high 
MeOH concentrations, the MeOH crossover and the mass trans- 
fer of air are mainly limiting the cell performance, especially 
at low and high current densities, respectively. Altogether, the 
developed model is in reasonable agreement with both the 
measured current distributions and polarization curves. The spa- 
tial information gained of mass transfer phenomena inside the 
DMFC is valuable for the optimization of the DMFC operating 
parameters. 
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